Accumulation and thermalization of cold atoms in a finite-depth magnetic trap 
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We experimentally and theoretically study the continuous accumulation of cold atoms from a 
magneto-optical trap (MOT) into a finite depth trap, consisting in a magnetic quadrupole trap 
dressed by a radiofrequency (RF) field. Chromium atoms (^^Cr) in a MOT are continuously optically 
pumped by the MOT lasers to metastable dark states. In presence of a RF field, the temperature 
of the metastable atoms that remain magnetically trapped can be as low as 25 fiK, with a density 
of lO^'^ atoms. m"'', resulting in an increase of the phase-space density, still limited to 7.10"^ by 
inelastic collisions. To investigate the thermalization issues in the truncated trap, we measure the 
free evaporation rate in the RF-truncated magnetic trap, and deduce the average elastic cross section 
for atoms in the metastable states, aei = 7.0 x 10~^^m^. We finally discuss the possibilities for 
using this scheme of continuous accumulation and RF evaporation to rapidly reach high phase-space 
densities from a MOT. 

PACS numbers: 32.80.Pj, 47.45. Ab, 32.80.Cy 
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Magneto-optical trapping is one of the greatest recent 
advances in atomic and molecular physics, opening many 
new areas in physics, including the study of Bose-Einstein 
condensation in dilute systems. One of the advantages of 
a magneto-optical trap (MOT), is that cooling and trap- 
ping are simultaneous. However, inherent limitations of 
the cooling mechanism, for example light-assisted colli- 
sions or multiple scattering of light, limit typical phase- 
space densities to 10~^ when strong fluorescence lines are 
involved in the cooling transition [ij. To reach higher 
phase-space densities, it is usually necessary to use se- 
quential operations, such as cooling in molasses, load- 
ing in conservative traps, and forced evaporation, which 
greatly limits the rate at which a Bose-Einstein conden- 
sates (BEG) can be produced. 

New possibilities arise for atoms whose electronic level 
structure includes metastable dark states (such as Sr, 
Er, Yb, or Cr). Spontaneous emission from the excited 
state of the cooling transition depumps atoms into a dark 
metastable state, which can be trapped, either optically, 
or magnetically. This provides an interesting means to 
continuously accumulate atoms in a magnetic trap 0] , or 
for a continuous loading of an optical or magnetic waveg- 
uide. Cold metastable atoms are directly produced inside 
the waveguide or the trap, which reduces heating associ- 
ated with the sequential loading from a MOT. This may 
have promising applications in the prospect of continu- 
ously producing a coherent beam of atoms, by perform- 
ing forced evaporation as the atoms propagate along a 
waveguide Q. 

In this paper, we study the accumulation of metastable 
chromium atoms in a magnetic trap (MT)jdressed by a 
radiofrequency (RF) field (see Fig [J) 0] (1 ■ As in 
and 0, chromium atoms in the '^83 ground state are 
captured in a MOT, and dcpumped to dark metastable 



states ^-D4,3. We observe different regimes as a function 
of the RF frequency. When the RF frequency is small, 
atoms in internal states adiabatically connected to high 
field seeking states at small magnetic fields are accumu- 
lated in a (3D) W-shaped potential. When the RF fre- 
quency i> is larger than the MOT temperature T^qt, 
i.e. hv > ksTMOT (where h is the Planck constant, and 
ks the Boltzmann constant), atoms in low-field seeking 
states are accumulated in a finite-depth potential. In this 
latter case, we achieve the continuous loading of a trap 
where evaporation is also continuous. 

The purpose of this paper is to carefully study the ac- 
cumulation of atoms in such a finite-depth trap, including 
the analysis of the different loss and heating mechanisms, 
the role of elastic collisions and the issues of evaporation 
and thermalization. We will show that fairly large phase- 
space densities can be reached. In the case of chromium, 
the gain in phase-space density is nevertheless limited by 
large inelastic loss rates associated with collisions with 
atoms from the MOT, and by large inelastic collision 
rates between metastable atoms. 

In the first part of this paper, we describe our exper- 
imental observations, which show the fast (300 ms) pro- 
duction of a cloud whose phase-space density can reach 
7.10"^. This is about a fifty fold increase compared to 
what is typically obtained in standard MOTs, and a ten 
fold increase in the phase-space density compared to the 
Cr MOT we operate. 

In the second part, we present a theoretical model 
based on rate equations for the atom number and the 
total energy of the sample, to interpret our experimental 
results and determine which processes limit the phase- 
space density. This model includes loading and loss pro- 
cesses, and assumes that thermal equilibrium is reached 
on a timescale which is short compared to the timescales 
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FIG. 1: Relevant level structure for chromium. The cooling transition 
is ^Ss — »^ P4. Atoms arc depumped to the metastable state. For 

this state, we also present a sketch of the RF dressed magnetic Zeeman 
eigen energies. Along the x axis, the avoided crossing occurs at RftF- 
The density profiles are also qualitatively represented (shaded area) for 
atoms in the MOT (a), and for metastable states depending on the sign 
of their magnetic quantum number (b) and (c). 

associated with loss processes. We also assume that the 
trap depth is large enough for the thermal distribution in 
the truncated trap to be identical to the thermal distri- 
bution in the non-truncated trap (Q, [^). The compar- 
ison between our experimental data and our simulation 
points towards a strong heating process due to collisions 
between the magnetically trapped atoms and the MOT 
atoms, which strongly limits the achieved phase-space 
density. 

In the third part, we describe additional experiments, 
which allowed us to measure the evaporation rate, by 
comparing the loss rates and the heating rates with 
and without RF. From these measurements, we infer 
the average cross section for elastic collisions between 
metastable chromium atoms, which was never measured 
before. This measurement indicates that thermalization 
is not achieved when the atoms are accumulated in the 
truncated magnetic trap at low RF frequencies. We fi- 
nally discuss the consequences of our observations in the 
prospect of reaching high phase-space densities by con- 
tinuously accumulating atoms in a truncated trap. 



CONTINUOUS ACCUMULATION IN A 
RF-DRESSED MAGNETIC TRAP 

Our experimental setup was described in a previous 
paper Chromium atoms are heated to 1500 C, de- 
celerated in a Zeeman slower, and captured in a MOT 
whose magnetic field gradient is 18 G/cm in the verti- 
cal direction, and 9 G/cm in the other two directions. 
Compared to reference , the MOT laser beam 1/e^ ra- 
dius is reduced to 3 mm. Here, we specifically study the 
accumulation of ^^Cr atoms in the state, in pres- 
ence of a RF field. We generate the RF field Bz using a 



commercial synthesizer, a 11 dB pre- amplifier, and a 10 
W amplifier, sending the RF power to a copper coil (2.5 
turns, 10 cm diameter). After a given time r {— 3s for 
most experiments reported in this paper) during which 
we run the MOT without any repumper, we shine the red 
repumping light (663 nm for ^D^ P3) for 4 ms, and 
capture an image of the cloud, using dark ground absorp- 
tion imaging along the horizontal y axis (1:1 imaging, 
using a telescope with two 2 inch diameter achromatic 
doublet lenses, and a CCD camera), with a 100 /is res- 
onant pulse well below optical saturation. Our imaging 
resolution is limited by the size of the camera pixels, 6.5 
fim. The repumping time, 4 ms, is long enough for all 
atoms in the state to be repumped to the ^83 elec- 
tronic ground state, but still short compared to all other 
timescales in the experiment. In addition, the recoil en- 
ergy associated with repumping is negligible compared to 
the kinetic energy of the atoms. We therefore consider 
that the temperature and densities measured for atoms 
in the ^83 after repumping are equal to those of atoms 
in the state before repumping. 

For in situ density measurements, the absorption pic- 
tures are taken in presence of the magnetic field gra- 
dient. We analyze the images by a fitting procedure 
which assumes thermal equilibrium in a trap potential 
V{x, y, z) = msgjfJ'BB' x'^ + + (m^ is the mag- 
netic sublevel, gj = 1.5 the Lande factor, /is the Bohr 
magneton, and B' = 9 G/cm the magnetic field gra- 
dient). From this analysis, we first deduce the total 
atom number. The cloud size, typically 500 /xm, is small 
enough that we can neglect the Zeeman effect on the cy- 
cling transition (the Zeeman shift is less than one tenth 
of the linewidth). 8ince the atomic sample is unpolarized 
and the magnetic field is inhomogeneous, we assume an 
average Clebsch Gordan coefficient of 3/7 for the reso- 
nant photon absorption cross section aabs — (where 
A = 425.5 nm is the photon wavelength), correspond- 
ing to an uniform magnetic sub-level distribution. From 
optical pumping calculations simulating an unpolarized 
sample, we estimate that the corresponding error bar on 
the number of atoms is ±15 percent: this is the main 
contribution to the error bar in the density and in the 
determination of chromium inelastic loss parameters de- 
scribed later in this paper. 

We measure the number of atoms accumulated in the 
truncated MT, as a function of time r, for different RF 
frequencies i'. At any frequency, we observe an exponen- 
tial loading of the trap, and we plot in Fig [5] the 1/e ac- 
cumulation time Tioad as a function of the RF frequency. 
The accumulation time increases approximately linearly 
from 0.3 s for ^ = 1 MHz to 1.4 s for > 6MHz. With- 
out RF, the accumulation time is 1.4 s. In Fig[3l we show 
the steady state atom number after 3 s of accumulation 
in the truncated MT as a function of the RF frequency. 

Since the atom cloud is not polarized, the ex- 
act analysis of the absorption images is complicated. 
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FIG. 2: Accumulation time in the truncated MT Tioad^ as a function 
of the RF frequency. 



For sake of simplicity, we will assume in this pa- 
per that all atoms experience an average potential 
msgj HbB' yaJ^+y^^^; where Ws is an average mag- 
netic quantum number, which we deduce from direct tem- 
perature measurements. The density profile at thermal 
equilibrium is umt = noea;p(— Vb/fcsTi/ + + 4z2), 
where T is the temperature, and Vq = msgjUBB' . From 
the analysis of the absorption images, we can deduce the 
1/e radius of the cloud d, as well as the peak atom density 
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, where O.D. is the experimental central op- 
tical depth. After 3 s of accumulation, the observed peak 
atom density is almost independent of the RF frequency, 
except when this frequency is set below approximately 1 
MHz (see FigE]). 
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FIG. 3: Atom number after 3 s of accumulation, as a function of the 
RF frequency (squares). Solid line : theory (sec text). 
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FIG. 4: Peak atom density after 3 s of accumulation (squares). Solid 
line : theory (see text). 



In addition to density measurements, we perform di- 
rect temperature measurement, by analyzing the free-fall 
expansion of the atoms after release from the magnetic 
trap. The current in the coils producing the magnetic 
field gradient is switched off in 500 ^is. Eddy currents 
generate a non-negligible bias magnetic field at the atom 
location for more than 5 ms, but this does not modify 
the ballistic expansion of the atoms. To deduce the tem- 
perature, we therefore analyze times of flight larger than 
600 /j-s, by fitting the integral of the absorption images 
by gaussian functions. Although the initial shape of the 
cloud is not gaussian, we numerically checked that the 
systematic error bar associated to this procedure is less 
than 10 percent, when fitting the vertical direction (along 
which the extension of the cloud is the smallest), for tem- 
peratures on the order of 60 /iK. 

Using this temperature and the measured 1/e size of 
the cloud in situ d, we deduce the mean magnetic quan- 
tum number ms ~ ksT / dgj^sB' of the atoms in the 
MT. Interestingly, does not change much with the 
RF frequency, staying close to 2. 

In figure [6l we show the steady state phase-space den- 
sity after r =3 s of accumulation, as a function of the 
RF frequency. Because the measured temperature of the 
atoms decreases with decreasing RF frequency (Fig [5]) , 
and because the measured peak atom density remains al- 
most constant (Fig|4|), the phase-space density increases 
at low RF frequencies (Fig [6]), to reach a maximum value 
of (7.0 ± 1.5) 10-6 at 1/ = 1.5 MHz. The error bar is 
due to the systematic errors on the density and on the 
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FIG. 5: Temperature of the atoms after 3s of accumulation, as a 
function of tlic RF frequency. Tlic experimental error bar on the tem- 
perature is about 10 /iK. Solid line : theory (see text). Experimental 
points at the left of the vertical dotted line correspond to a situation 
where thermal equilibrium is not reached (sec text). 



temperature measurements (see above). The maximum 
phase-space density value is significantly higher than typ- 
ical phase-space densities achieved in a MOT involving a 
strong resonance line. 

For frequencies lower than 1.5 MHz, this trend inverts: 
when lowering the RF frequency, the density of the atoms 
decreases (see Fig|4|), and we infer a decrease of the phase- 
space density, instead of an increase. For such low RF 
frequencies, the atoms in the low-field seeking states see 
a trap whose depth is lower than the MOT temperature, 
and most of them are therefore not trapped. However, for 
frequencies on the order of 800 kHz or below, new inter- 
esting features appear in the absorption images. The 1/e 
size of the imaged cloud becomes much larger, and the 
atom lifetime is longer than the lifetime measured at 1 
MHz. We interpret these pictures as resulting from atoms 
trapped in the W-shaped adiabatic potential correspond- 
ing to states, dressed by the RF field, adiabatically con- 
nected to high- field-seeking states at B=0 Indeed, 
the cloud size increases linearly with RF frequency, and 
its vertical position also shifts linearly with RF frequency, 
as is expected in the RF-dressed potential, whose size is 
determined by Rhf = hi^Rp/gjf^BB' (see Fig. ([T])). For 
low RF frequencies, we found that the number of atoms 
increases before decreasing for large RF power. We at- 
tribute this effect to higher harmonic frequencies appear- 
ing at high RF power : a substantial power of the second 
harmonics of the RF frequency is then exciting the atom 
sample, so that the atoms reaching 2Riip have a large 
probability of being extracted from the dressed trap. 

In the remaining of this paper, we will focus on RF fre- 
quencies larger than 1 MHz, for which we only trap low- 
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FIG. 6: Phase-space densities after 3 s of accumulation in the RF 
trap, as a function of the RF frequency (Squares). Solid line: theory. 
Experimental points at the left of the vertical dotted line correspond 
to a situation where thermal equilibrium is not reached (sec text). 

field seeking states. One of the main goals of this paper 
is to quantitatively understand the achieved phase-space 
densities in this regime. In the next part, we describe a 
theoretical model, which assumes that the elastic scat- 
tering rate is high enough to insure thermal equilibrium. 
Then, in Part 3, we will experimentally determine the 
elastic scattering rate, in order to determine in which 
RF regime this theoretical model is expected to be valid. 



THEORETICAL MODEL AND 
INTERPRETATION 

Our theoretical model for the accumulation of the 
atoms in the RF-truncated magnetic trap is based on 
rate equations for the number of atoms and the total 
energy of the system. It is elaborated from Q and [§]. 
We include in our model the density dependent loss pro- 
cesses, assuming that the density profile of the atoms is, 
at all time, well described by the Boltzmann distribu- 
tion at thermal equilibrium in the non-truncated trap. 
This relies on two main assumptions: first, the elastic 
collision rate is large enough to reach thermal equilib- 
rium on the timescale of the experiment (defined by the 
dominant loss process); second, the trap depth is large 
enough compared to the magnetic trap temperature, so 
that one can neglect modifications in the thermal distri- 
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bution linked to the finite depth of the potential (01, [I])- 
In our experiment, thermal equilibrium is most likely not 
reached at the beginning of the loading of the MT, when 
the number of atoms is still small. For this reason, we 
do not expect this model to describe the dynamics of the 
loading. However, the steady state reached at the end 
of the accumulation is the one predicted by our model, 
provided the rate of elastic collisions in the state is 
high enough (or if the steady state density profile is close 
enough to the Boltzmann density distribution). 

To be able to capture the physics of the accumula- 
tion in the truncated MT, our model takes into account 
the loading of the atoms and various loss mechanisms. 
The key point for the loading parameters is that the 
loading rate and the average energy per atom loaded in 
the MT both depend on the RF frequency. We exper- 
imentally checked that the MOT atom number is not 
modified by the presence of the RF field. Wc there- 
fore assume that the production rate of atoms in the 
metastable states does not depend on the RF frequency 
V. However, once the atoms decay into the ''^D4 state, 
we assume that they are only trapped if their total en- 
ergy is smaller than the trap depth. The loading rate 
T{v) is therefore obtained by multiplying F, the loading 
rate without RF, by the probability P{v) that an atom 
transferred from the MOT - at a position (a;,j/, z), and 
with a velocity v - into the metastable state has a total 
energy Etot = l/2mw^ -I- V^yJ x'^ + y'^ + Az'^ smaller than 
the depth of the truncated MT, rnshv: 



(1) 



P(i') can be evaluated numerically by integrating the 
phase-space density over the domain D defined by the 
condition Etot < rnshv : 
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Our model includes different loss mechanisms: 

(i) collisions with the background gas, mostly hot 
atoms coming from the oven (rate Fq); 

(ii) losses corresponding to inelastic collisions with 
atoms from the MOT, with a rate 



Fi = 1/iV j 



(4) 



where N is the total number of atoms in the magnetic 
trap, and Ppjj is the inelastic loss parameter for collisions 
between D atoms and P atoms; 

(iii) two-body losses (rate F2 = L2no/2^, where uq is 
the MT peak atom density); 

(iv) three-body losses (rate F3 = Lsnl/S^. In general 
L3 has to be measured. However, in the limit of the 
existence of a weakly bound dimer, L3 is related to the 
scattering length of the atoms to the fourth power [ld|); 

(v) Majorana losses (rate Tmaj — C^'aJII^^ where /i 
is the atom magnetic moment, A/i is the change in mag- 
netic moment due to the Majorana spin flip, C is a nu- 
merical factor, on the order of 3, and d is the 1/e radius 



of the cloud [11|) 



(vi) Evaporation related to elastic collisions, with rate 
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In this formula, we have used the fact that the 
phase-space density is proportional to the spatial den- 
sity in the MOT, nMOT(^) = "-o,mot exp(-x^ + y'^ + 
Az'^)/2w\jQrj,), and to the -y/uexp(— it) factor coming 
from the free space kinetic energy distribution. This 
MOT density mostly reproduces the experimental one, 
and allows great mathematical simplifications. We found 
in our simulations that the exact form of the MOT den- 
sity mainly affects the loading parameters and the results 
of our simulations at low RF frequencies, but doesn't 
modify much the inelastic loss parameters that we de- 
duce from the simulations. 

We also calculate the average total energy per atom 
trapped in the truncated MT, E{v): 



where Fg; — noUeiv is the elastic collision rate, with 
a el the elastic cross section and v — \ ^^-'-^^ the mean 

y 7rm 

thermal velocity. 77 is the ratio of the trap depth to the 
trapped atom temperature, and /(??) is given in equation 
(11) of ref ^81] (which is valid for ?7 > 4). 

We neglect any elastic collisions between atoms in the 
metastable states and atoms from the MOT 1^. The 
atom number and total energy rate equations read: 



dN 
~dt 



= T{v) - (Fo + Fi + F2 + F3 + F„,aj + Fe„)iV (6) 



d{%/2NT) 
Jt 



= E{v)T{v)- /'^'^^ (^) 

2— 0,l,2,3,maj,ei) 



where the fikpT's describe the average energy loss 
when one atom is lost by the corresponding mechanism. 
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In a linear trap, one shows that /o = 9/2, /2 = 3, 
/a = 5/2, fmaj = 5/2. For instance, to evaluate /2, 
one needs to integrate over the whole cloud the local in- 
elastic rate L2"mt ('') times the local mean energy of the 
atoms V{r) + 3/2kBT: 

f2T2NT = -L2 J dVn2,^(r) {V{r) + 3/2kBT) (8) 

If the trap is deep enough, f^v — rj + 1 (see In 
addition, 

fiTiNT = j (frPpDnMOT{r)nMT{r){V{r) + 'd/2kBT){9) 

We stress that there are few adjustable parameters, 
since uq^aiot = 7.5 x 10"^^ m~^, Tmot = 120 /uK and 
wmot — 100 ^m are independently measured. As shown 
below, we determine L2 from an independent measure- 
ment. T3 and Tev both depend on the scattering length 
a, which is used as a free parameter here. F and l3pD are 
the two other free parameters. 

We numerically solve equations (O and ([7]) using the 
experimentally measured ms for each RF frequency. We 
first fit our numerical results to the experimental results 
for large RF frequencies {v > A MHz). In this regime, we 
expect that evaporation is negligible, because the trap 
depth is much larger than the temperature of the atoms: 
for instance, a RF frequency of 8 MHz corresponds to 
•q > 10, so that the evaporation time is much longer than 
the accumulation time in the truncated trap. Therefore, 
we set Fet, = for the comparison to the experiment at 
large RF frequencies. The comparison to our experimen- 
tal results for the temperature and number of atoms after 
3 s of accumulation set the parameters F = 3.3 x 10^ 
and (ipD = (4.9 ± 0.3 ± 0.5) x 10-^° cm^s-i. The error 
bars are respectively the statistic error bar of the fit, and 
the systematic error bar associated to the uncertainty on 
the total atom number. With these values, the agree- 
ment between theory and experiment is satisfactory, for 
RF frequencies exceeding 4 MHz (see Fig [3l SI [5] and [6]) . 
We note that the value oi (ipo that we deduce from our 
analysis is close to what was measured in reference [isj . 

For large RF frequencies, we checked that evapora- 
tion, Majorana losses, and three-body losses are negli- 
gible. The only other parameter which is not negligible 
is the inelastic two-body loss parameter, which is mea- 
sured independently in our experiment (see below). The 
main conclusion of our analysis at high RF frequencies is 
that inelastic collisions with atoms from the MOT limit 
the total number of atoms that are accumulated in the 
MT, and produce a strong heating: the temperature of 
the atoms in the MT, which, in absence of collisions, 
could be as low at Ti\iot/S [l], raises to about 100 /xK for 
large RF frequencies. This drastically limits the achieved 
steady state phase-space density. 



At RF frequencies lower than 4 MHz, the role of evap- 
oration is not negligible. However, we find that varying 
the elastic cross section aei does not lead to large modi- 
fications of the numerical results, so that it is in practice 
impossible to deduce from a fit a precise value of the elas- 
tic scattering rate. In addition, for the lowest frequencies, 
our model, which relies on the assumption that the gas 
is thermalized to the Boltzmann distribution, may not 
be fully applicable. The loading time Tioad of the atoms 
in the truncated MT gets smaller and smaller as the RF 
frequency decreases (see Fig. [5]): at low RF frequencies, 
the spatial overlap between the MT and the MOT gets 
better and better, which increases the loss rate related 
to inelastic collisions with atoms from the MOT (process 
(ii)), and reduces Tioad- Therefore, Boltzmann thermal 
equilibrium is not necessarily reached before the steady 
state. This will be studied in detail in the following sec- 
tion. 
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FIG. 7: ?7 as a function of the RF frequency. For 17 < 4, the density 
distribution at thermal equilibrium is modified by the trap truncation, 
which is not taken into account in the theoretical model. 



Effects linked to the small trap depth (and their con- 
sequences on the thermal distribution) may also have an 
important role at low RF frequencies. Our model, assum- 
ing that the thermal distribution in the truncated trap 
is identical to the one in a infinite depth potential, will 
only be valid provided 77 > 4 Q. Although the sample is 
not polarized, the measurement of the 1/e radius of the 
cloud d provides a good estimate for 77, the ratio of the 
trap depth [mshv) to the trap temperature ksT = V^d, 
because this ratio does not depend on ms- We plot 77 in 
Fig [71 and find that our theoretical model may not be 
applicable when the RF frequency is smaller than about 
2 MHz. 
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EVAPORATIVE COOLING AND 
THERMALIZATION IN A RF-TRUNCATED 
MAGNETIC TRAP 



To investigate thermalization issues, we perform a new 
set of experiments to measure the lifetime and the heat- 
ing rate of atoms in the ^1)4 state with and without 
the RF field. For these experiments, the atoms are first 
accumulated without RF for 3 s in the magnetic trap. 
Then, the MOT is switched off, and the number of atoms 
and temperature of the cloud are measured after an ad- 
justable delay t. This second step of the experiment is 
first realized without RF, then in presence of a 3 MHz 
RF field. 

When no RF field is applied, we measure the decay 
of the metastable atoms in the MT (see Fig [5]), from 
which we deduce the inelastic loss parameter L2 (see Fig 
[51 which shows a clear non-exponential decay) . We fitted 
the data of Fig [5] assuming, in addition to a one-body 
loss, a two-body loss parameter. Given the low atomic 
density, we assume that three-body processes are neg- 
ligible. During the decay, the cloud volume increases 
approximately linearly. In the fit of Fig [51 we assume a 
linear dependence to the volume on time, based on a fit 
to the measured time-dependent size. We found the in- 
elastic loss parameter to be Lo = (3.3±0.5±0.5) x 10~^^ 
cm^s~^ (in agreement with [13[). The respective error 
bars are statistic error bar of the fit, and systematic error, 
mostly coming from the uncertainty on the total number 
of atoms. We stress that this inelastic loss parameter is 
much larger than the expected dipolar relaxation rate, 
which rules out dipolar relaxation as the main inelastic 
loss channel 
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We repeat this lifetime experiment in presence of a 3 
MHz RF field during the decay. Then, the decay of the 
cloud is modified by evaporation, which changes both the 
total loss rate, and the evolution of the temperature (see 
Figinj). Our goal is to use these observations to estimate 
the evaporation rate and infer the elastic collision rate. 
However, one must first evaluate the modification of the 
lifetime of the cloud linked to spilling of atoms out of the 
trapping volume set by the RF frequency: the experi- 
mentally observed heating produces atoms whose energy 
is larger than the trap depth. They then leave the trap, 
without undergoing elastic collision, assuming that the 
mean free path of the atoms is larger than the size of 
the sample. Heating without RF translates into spilling 
when RF is applied, and the rate at which spilling oc- 
curs is therefore linked to the heating rate without RF. 
We measure the time evolution of the I/e radius of the 
cloud d without RF to estimate Tspui ~ 0.1 s~^, the 
rate of atoms leaving the trapping volume set by the RF 
frequency, i.e. the rate of atoms spilled out of the trap 
volume. 

From the atom number decay at short times, we de- 
duce the experimental values of the atom loss rates with- 
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FIG. 8: Decay of atoms in the D state (squares). Solid line: result 
of the fit (sec text). The dotted line is a result of a purely exponential 
fit of the experimental data for t> 5s. The dashed line is a result of an 
exponential fit for t< Is. 



out RF (see Fig [8]), Tioss = l/T,oss = 1/2.4 s-\ and 
with RF, Tioss,RF — 1/1.35 s^^. We can therefore esti- 
mate the value of Fei, -I- TspUi = Tioss,rf - Tioss = 0.32 
s^^ . From this equality, we obtain a first estimate for the 
evaporation rate from the modification in the lifetime of 
the cloud in presence of RF: T^v ~ 0.22 s~^. 

We deduce an additional independent estimate for the 
evaporation rate by comparing the heating rate in the RF 
truncated trap to the heating rate without RF. With the 
3 MHz RF field on, we observe no substantial heating, as 
proven by temperature measurements as the hold time t 
is varied (see Fig[n|). This indicates that the cooling rate 
due to evaporation balances the heating rate related to 
two-body inelastic collisions. For RF frequencies higher 
than 3 MHz, the cloud does heat, whereas we observe 
cooling when the RF frequency is lower. This observa- 
tion allows us to make a second determination of Fe^. 
Using ([7|), assuming only one-body, two-body losses and 
evaporation, and taking into account the experimental 
observation that ^ ~ 0, we have: 



9 dN 

^-^T « -/oiYTFo - hNTT2 - {v + l)TevNT (10) 



from which we deduce, using eq.®: 



(9/2-/2)r2 



r? + 1-9/2 



(11) 
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FIG. 9: Time evolution of the temperature without RF (squares) and 
with a ly — 3 MHz RF field (triangles). Lines arc guides for the eye. 
The error bar illustrate the typical combined statistical and systematic 
errors on temperature measurements. 



We estimate r2,t=o = Tioss - Tq, where Tq = I/Tq = 
1/9 is the loss rate at large time delays t, for which 
the atom number is smah enough to neglect inelastic col- 
lisions (see Fig ^ . We attribute Fq to the effect of col- 
lisions with the background gas (mostly hot atoms com- 
ing from the oven). The 1/e size of the cloud and the 
experimentally measured temperature give ms ~ 3.5, 
and an average trap depth to trap temperature ratio of 
rj w 5.5 ± 1. We therefore deduce from eq. pTjl that 
Fe„ = 0, 2 ± 0.04 s~^, in reasonable agreement with our 
first estimate. 

From these measurements, we deduce the elastic colli- 
sion rate, using eq. ([5]) and eq. (11) of ref [sj to estimate 
the evaporation fraction f{ri): T^i = 20 ± 4 ± 11 s""'^ (er- 
ror bars respectively correspond to statistical error on Fe„ 
and to the estimated systematic error on rj). The experi- 
mentally measured temperature is 100 iiK, and the peak 
density 10^^ cm~^, at t = 0. We therefore infer an aver- 
age elastic cross section of =(7.0 ± 1.4 ± 3.5 ) x 10"^^ 
m^, close to the unitary limit at this temperature. The 
error bars respectively correspond to the statistical er- 
ror, and the estimated systematic error (dominated by 
the estimate of 77). This is, to our knowledge, the first 
measurement of an elastic cross section for a transition 
element in a metastable state. It should be mentioned 
that, in relation with the prospect for Bose-Einstein con- 
densation with alkaline earth elements, the elastic cross 
section was also recently determined in metastable Ca 

Using this average collision rate, we now check if our 



theoretical model can be used to analyze the evaporation 
during the decay, i.e. if thermal equilibrium is achieved 
throughout the decay, so that we can use equation pT|) . 
To be able to describe the cloud by the Boltzmann dis- 
tribution throughout the decay, it is necessary that a 
sufficient number of elastic collisions take place during 
the typical evolution time of the cloud. It is difficult to 
define exactly a timescale to reach thermal equilibrium, 
as relaxation to the Boltzmann distribution is not expo- 
nential when this relaxation is driven by intraparticles 
collisions (as can be seen for example from eq. (14) in 
Q). We perform numerical simulations of the Boltzmann 
equation, and find that, in practice, a ratio of about 4 
elastic collisions per inelastic collision is enough to main- 
tain the energy distribution of the most populated states 
close to a Boltzmann distribution. In our experiment, 
66. The distribution should therefore be 
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close to the Boltzmann distribution, and our determi- 
nation of the elastic scattering rate in the ^Z?4 state is 
therefore self-consistent. 

We now can turn back to the analysis of the experi- 
ments described in the first part of this paper, and to 
the issue of thermal equilibrium in these experiments. 
For the continuous accumulation in the RF truncated 
trap, the loading time is mostly set by inelastic collisions 
with the MOT atoms. For low RF frequencies, for which 
the spatial overlap between the MT and the MOT is the 
best, the loading time Tioad is short, between 300 ms and 
Is (see Fig[2|). Therefore, the elastic scattering rate is 
too small to provide thermal equilibrium at the end of 
the loading process. In Fig 1101 we plot the approximate 
number of collisions before the steady state of accumu- 
lation is reached, qualitatively given by UMTO'eivTioad/'^- 
From Fig. [TOl we expect a lack of thermalization at RF 
frequencies below 4 MHz. 
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FIG. 10: Number of elastic collision before the steady state is reached 
when atoms are accumulated in the RF-truncated magnetic trap. When 
this number is smaller than approximately 4, thermal equilibrium is not 
reached, and our theoretical model does not apply. 



Preliminary numerical simulations based on the Boltz- 
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mann equation show that indeed thermalization is not 
achieved at low RF frequencies. For higher RF frequen- 
cies, these simulations also show that thermalization is 
not perfect. Since thermalization is not achieved in the 
experiment, especially at low RF frequency, the agree- 
ment between experimental results and theory may look 
surprisingly good. Such a good agreement relies on the 
fact that, at each RF frequency, the experimental den- 
sity profile is always well fitted by a density profile assum- 
ing thermal equilibrium in a potential Vq ^/ x'^ + y"^ + 4z^. 
Therefore all loss terms in the system are well described 
by the terms Fj and fiTi in equations ^ and ([7]) , which 
therefore still describe the experimental system, even 
if thermal equilibrium is not reached. For these rea- 
sons, the measurements of the inelastic loss parameters of 
chromium given in this paper should be correct. Further- 
more, even if technically temperature and phase-space 
densities are not defined at the left of the vertical dot- 
ted line in Fig [5] and \6\ the quantities we plotted still 
have a relevant physical meaning: we measured the spa- 
tial densities and the total energy of the system, which 
will define the thermodynamics quantities of the sample, 
once thermal equilibrium is reached. 

In this paper, we have studied the accumulation of 
metastable atoms in a magnetic trap of finite depth. We 
demonstrated that phase-space densities 50 times larger 
than the typical ones in a MOT can be reached in less 
than Is. However, we showed that inelastic collisions 
with atoms from the MOT reduce the loading time of 
the atoms in the trap so much that thermalization is not 
reached when the trap depth is too small. In addition, 
the heating rate associated with these inelastic collisions 
greatly reduce the achievable phase-space densities. This 
feature may be specific to chromium. With atoms for 
which such a loss process is reduced, we suggest that the 
continuous accumulation of atoms in a finite depth trap 
(either a magnetic trap or a dipole trap) could be an 
efficient way to rapidly reach high phase-space densities. 
Finally, we also studied the inelastic collision properties 



of metastable chromium atoms, as well as their elastic 
properties, providing the first measurement of the elastic 
scattering rate for chromium atoms in the state. The 
ratio of the elastic rate to the inelastic rate is close to 
five, ruling out traditional forced evaporation techniques 
to reach BEG in these states. 
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